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A small number of laser-cooled ions trapped in a linear radiofrequency multipole trap forms a 
hollow tube structure. We have studied, by means of molecular dynamics simulations, the structural 
transition from a double ring to a single ring of ions. We show that the single-ring configuration has 
the advantage to inhibit the thermal transfer from the rf-excited radial components of the motion 
to the axial component, allowing to reach the Doppler limit temperature along the direction of the 
trap axis. Once cooled in this particular configuration, the ions experience an angular dependency 
of the confinement if the local adiabaticity parameter exceeds the empirical limit. Bunching of the 
ion structures can then be observed and an analytic expression is proposed to take into account for 
this behaviour. 



I. INTRODUCTION 

In the wide range of radiofrequency (rf) traps, 
quadrupole traps take a particular place due to the 
parabolic rf potential they create and to the harmonic 
pseudo-potential it generates. In contrast, higher order 
traps show a more "rectangle" shape of the rf electric 
field which induces a flat pseudo-potential well with a 
reduced rf-driven micromotion in the center of the trap 
[l[ . This feature finds its application in the spectroscopy 
of cold molecular samples |2| or in the interrogation of 
the hyperfine transition of a large atomic sample for fre- 
quency metrology 0] , the ions being cooled by collisions 
with a buffer gas. 

Laser-cooling of ions trapped in multipole traps has 
only been performed very recently Q and experimen- 
tal observations have demonstrated a tube-like structure 
for the ion organization Q. This self-organization in 
empty core structure results from the competition be- 
tween Coulomb repulsion and a confining force very soft 
in the center of the trap and stiffer towards the walls. 
For large sample, this behavior is confirmed by the den- 
sity profile deduced in the cold charged fiuid limit [Qj and 
for smaller ones, by simulations of the equilibrium struc- 
tures 0, Si • For a large number of ions trapped in a 
linear multipole trap, the general structure of a crystal- 
lized sample is a set of concentric tubes for which the 
revolution axis is the longitudinal trap axis. The evo- 
lution of the structure with the number of ions and the 
trapping parameters is extensively studied in @ . 

As an example, figure [1] shows the simulated structure 
of 1000 cold calcium ions inside the pseudo-potential of 
a linear octupolc trap, laser-cooled to the Doppler limit 
temperature of 0.5 mK. For the same trapping parame- 
ters and a smaller number of ions, the structure reduces 
to a single tube, formed of parallel rings of ions aligned 
along the trap axis. For a decreasing number of ions, the 
tube reduces to two rings then to a single ring. Ring- or 
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FIG. 1. (color online) Structure of a set of 1000 laser-cooled 
ions in the pseudo-potential of a linear octupole trap. The 
structure is composed of three concentric tubes extending 
along the trap axis, which is in the z direction. The tem- 
perature is 0.5 mK in the three directions of space. 



tube-like structures offer the possibility to study dynam- 
ics in a two dimensional configuration. Furthermore, ion 
rings have symmetry properties that can be exploited in 
high-resolution spectroscopy. Among these highly sym- 
metric structures, the single ring configuration has the 
outstanding property to be a true two dimensional sys- 
tem where the motional temperature along the symmetry 
axis can reach the Doppler limit when the ions are laser- 
cooled, despite the rf-driven motion in the transversal 
plane Q. A previous numerical study [l^l of the cooling 
properties of ion structures formed in quadrupole traps 
has shown that for structures of reduced spatial dimen- 
sions. Coulomb interaction does not always couple the 
different directions of motion of the ions. This decou- 
pling inhibits the thermal transfer between the decoupled 
directions of the motion and allows to reach the Lamb- 
Dicke regime for the ion motion along the longitudinal 
trap axis @- 

In the present article we use molecular dynamics sim- 
ulation of laser-cooled ions in a linear octupole trap to 
study their structural phase transition from a two-ring 
to a one-ring configuration. The model used to simu- 
late the ion dynamics is described in section |TI] together 
with some brief reminders on linear multipole traps. Sec- 
tion mil presents an analysis of the structural transition 
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itself and enlightens the role of the rf-driven motion on 
this transition. In section IIV) we extend our study by 
increasing the adiabaticity parameter rjad {r) beyond the 
empirical limit. 



II. FRAMEWORK OF THE MODEL 

The numerical simulations use molecular dynamics, 
based on the symplectic velocity- Verlet algorithm to in- 
tegrate Newton's equations. We have simulated the dy- 
namics of a set of 100 ions in the potential of a linear 
octupole trap, cooled by Doppler laser-cooling. Ions are 
thus subject to conservatives forces deriving from the 
trap electric potential and the Coulomb inter-ion poten- 
tial, and dissipative forces due to Doppler cooling. The 
general form for the ideal rf potential of a linear 2fc-pole 
trap can be written in polar coordinates like 



$2fc ir,ip,t) = qeV2k — cos (kip) cos {ilt) , (1) 

where V2k is the amplitude of the rf potential applied to 
the 2k rods, oscillating at frequency r2/27r, tq is the in- 
ner radius of the trap, and qe is the electric charge of the 
ions. Unless otherwise specified, we use the potential of 
an octupole trap {2k = 8) to run the simulations. This 
rf potential induces an rf-driven motion which produces 
a ponderomotive force, which leads to the radial con- 
finement of the ions. It can be derived from the static 
pseudo-potential [l|: 
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The motion in this pseudo-potential has a time scale 
larger than an rf period and is called macro-motion. 

The axial confinement is achieved by a harmonic static 
potential well characterized by the frequency u!z/2Tr, 
which scales like -y/ qe/m. Obeying Laplace's law, this 
static potential has a deconfining contribution in the ra- 
dial plane. Hence, for each ion of mass m and position 
(r, z), the total trapping potential Utrp writes: 
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When the order of the multipole is higher than 2k = 
4, the deconfining contribution of the axial confinement 
leads to a shift of the minimum of the trapping potential 
in the radial plane, from the centre to a circle of radius 
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Hence, contrary to a quadrupolc trap configuration, the 
local minimum of the trapping potential is no more super- 
posed to a node of the rf electric field Throughout 



this article, the confinement parameters and number of 
ions are such that the ions sit very close to the potential 
minimum and they do undergo micro-motion. 

The Coulomb interaction is taken into account by the 
total pair interaction potential, which writes: 
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where Np is the number of ions and Sij 



is the relative dis- 
tance between ions i and j. The Doppler laser-cooling is 
modeled using velocity kicks simulating the ion's recoil, 
which happens depending on the probability laws of ab- 
sorption (given by the stationary solution of the optical 
Bloch equations) and emission (given by the exponential 
decay law). For symmetry reasons, laser-cooling is ap- 
plied in the three directions of the trap. Details about 
the model we used for the cooling are given in [lo| . 

Initial conditions are prepared by giving to each ion an 
initial position randomly chosen in the inner 10% of the 
radius of the trap and a random velocity obeying a ther- 
mal distribution corresponding to 1 K. Ions are trapped 
by the rf electric field and are subject to Doppler cool- 
ing. The system is ready for further simulations when 
the temperatures are stationary. The temperature is one 
of the relevant parameters to get insight into ion dynam- 
ics. As the rf-driven motion is not a thermal motion, it 
has to be subtracted from the velocities before comput- 
ing the temperature [H, In practice, this is done 
by averaging the velocity of the ions over an rf period. 
Each averaged component of the velocity is then taken 
into account to define the temperature associated to each 
direction of the motion. Explicit details about the com- 
putation of the temperatures can be found in [lo| . 

All the forthcoming results assume perfect trap oper- 
ation and neglect defects of the potential due to patch 
potential and/or misalignments of the electrodes as well 
as collisions with the residual gas. Patch potentials with 
a highly non-uniform spatial repartition could destroy 
the highly symmetric structure formed by few ions in a 
linear multipole trap. For high precision spectroscopy 
in multipole trap, more extra electrodes than for a line 
in quadrupole trap may be required to compensate for 
patch potentials and their gradient. Misalignment of the 
electrodes has been identified in a 22-pole trap IJ] as 
responsible for 10 extra potential minima, which cause 
bunching of the ions in these local potential wells. As 
for collisions, they may perturb the formation process of 
the crystal but the trapping conditions chosen for these 
studies imply a strong confinement in the axial and in the 
radial direction, of the same order of magnitude as the 
one met in quadrupole trap. Therefore, the stationary 
structures described in this article should survive to col- 
lisions like ion chains do. This strong confinement may 
seem unrealistic regarding the smooth profile of the trap- 
ping field in the center of the trap. The scaling law for 
the rf electric field in an octupole is Vs/^o^, and to reach 
a Rmin value of 120/xm with = 1 MHz and = 20 
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MHz like chosen in the foUowing, it takes an amphtude 
Vs equal to 208 V in a 0.5 mm inner radius trap, which 
is very accessible. As for quadrupole trap dedicated to 
single ion experiment, a multipole trap able to keep a 
single ring stable has to be specifically designed to that 
purpose. 



III. DOUBLE-RING TO SINGLE-RING 
STRUCTURAL TRANSITION 

The trapping parameters are such that the Coulomb 
repulsion is negligible compared to the deconfining ef- 
fect associated to the axial confinement and the radial 
position of the ions is controlled by the potential min- 
imum i?inin. Depending on the number of ions Np and 
the applied axial potential, the stationary structure is a 
single or a double ring. By stability consideration, one 
can show that there is a critical radius i?crit below which 
the equilibrium configuration is a double ring Q. This 
critical radius only depends on Np and uj^ through 
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For a set of 100 ions and uj^jl-K = 0.5 MHz, the stability 
critical radius i?crit is 181 /im. The rf electric potentiel is 
chosen such that the radius i?min is 120 /im and the ions 
organize themselves in two rings of 50 ions separated by 
an axial distance of about 9 /im. This double-ring struc- 
ture can be seen on figure [51 for an rf phase correspond- 
ing to the maximal amplitude of the micro-motion. The 
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FIG. 2. (color online) Double-ring structure of 100 laser 
cooled ions, in the rf potential of an octupole trap. The ax- 
ial confinement frequency is a;z/27r = 0.5 MHz, leading to an 
axial separation between the two rings of about 9 /im. The 
picture corresponds to the instant of the maximal amplitude 
of the micro-motion. 

pattern formed by the two ion rings is equivalent to an 
annular zig-zag structure [isj , the positions of two con- 
secutive ions being on different sides of the z = plan. 
This perfect annular zig-zag is obtained by choosing an 
even number of ions. In this configuration, the tempera- 
ture associated with the axial component of the motion is 
Tz~\ mK (sec figure [3]). If the velocity of the ions is ex- 
pressed in cartesian coordinates, the and Ty tempera- 
tures are identical and arc approximately equal to 12 mK. 



The difference between these values is attributed to a par- 
tial thermalization of the degrees of motion, the higher 
temperature being induced by rf-driven motion in the ra- 
dial plane. This motion induces rf-heating by which the 
kinetic energy of the rf-driven motion is transferred to the 
macro-motion [l^. It also reduces the Doppler cooling 
efficiency by the Doppler shift on the cooling transition. 
The temperature discrepancy shows that the Coulomb 
coupling in the two-ring structure is not strong enough 
to thermalize the motion in the radial plane with the 
motion along the axis [l^l . Inside the radial plane itself, 
the use of cylindrical coordinates reveals an important 
discrepancy in the temperature associated to the radial 
motions (T^ ~ 10 mK) and the one associated to the 
tangential motion (Tj^ ~ 2 mK, see figure [3|). 
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(color online) Temperatures associated to the mo- 
a polar coordinate system T,- (black upper curve), 
(red medium curve) and (green bottom curve), when 
the axial confinement frequency cj^ is increased. The vertical 
dashed line marks the theoretical transition frequency from 
the double-ring to the single-ring structure. The horizontal 
blue line materializes the Doppler limit temperature. 

Starting from this configuration, we study the double 
to single ring transition. To that purpose, the character- 
istic frequency of the axial confinement lo^I^t:^ is slowly 
increased from 0.5 MHz to 1.5 MHz. At the same time, 
the amplitude of the rf electric field is raised to balance 
the enhancement of the static deconfining effect in the 
radial plane, keeping the minimum of the trapping po- 
tential, i?min, and thus the radius of the ring, constant 
(see equation [4]). For a constant number of ions, the 
structural transition is expected from equation [6] to take 
place for uj^jl-K = 0.93 MHz. It is verified that the vari- 
ation of the trapping parameters is sufficiently slow to 
consider the system to be in a stationary state during the 
simulation. Particular cares are taken for values above 
0.92 MHz for which the increase of the confinement pa- 
rameters is halted every 0.1 MHz, by time periods of 
10 ms. 

The structural transition is first monitored by the dis- 
tance between the center of mass (cm.) of the two rings, 
Azcm, given on the left axis of figure HI as a function of 
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the characteristic frequency uoz/iir. The value for which 
the two rings merge is in good agreement with the the- 
oretical frequency of 0.93 MHz, for which the structural 
transition is predicted. 
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FIG. 4. (color online) Relative distance between the cm. 
of the two rings, Azcm (red curve, left axis), and standard 
deviation of the axial component of the ion's positions from 
the cm. position of the ring of which they belongs, Uz (black 
curve, right axis), function of the increasing axial confinement 
frequency ojz. The vertical dashed line marks the theoretical 
transition from the double-ring to single-ring configuration. 

The standard deviation of the axial component of the 
ions' positions from the cm. position is also plotted on 
figure 131 For uj^/^tt = 0.5 MHz, the standard deviation 
Uz is due to the axial thermal motion around the equi- 
librium position. Increasing lOz, the axial temperature 
(shown on figure [3]) drops while the local steepness of the 
potential in the axial direction is increased. Hence, the 
axial dispersion on ions' positions which is due to ther- 
mal motion is expected to decrease when ujz grows. Fig- 
ure |4] shows that it becomes maximum around the struc- 
tural transition before decreasing to a lower value than 
the initial one. To understand the transition process, 
we focus on the ion organization just before the struc- 
tural transition. While the two rings are getting closer, 
a distortion of the ring structures progressively appears, 
visible on figure [5] and stationary over the rf period. The 
distortion does not affect the equilibrium radius of the 
structures but the axial inter-ion distance is smaller for 
ions localized close to the rf electrodes than for those 
which are trapped far from the electrodes (with the con- 
vention defined by equation [Tl electrodes are at angular 
positions <p„ = n7r/4). This behaviour is attributed to 
micro-motion whose amplitude e)j (r, (ys, t) has the ana- 
lytic form: 

e), (r, t) = — — (cos (4(^) f - sin (4(^) ^) cos (Vlt) , 

(7) 

When ions are localized in front of an rf-electrode, the 
micro-motion amplitude is almost exclusively radial and 
has no significant consequence on the relative distance of 
two neighbouring ions belonging to different rings. On 



the contrary, when ions are localized at an inter-electrode 
angular position, the micro-motion is mostly tangential 
to the ring structure and depends like sin(4</j) on the 
ion angular position. Two neighbouring ions, separated 
by an angle 2Tr R/Np, do not have the same tangential 
amplitude of micro-motion, producing an oscillation of 
their angular separation. Since the neighbouring ions do 
not belong to the same ring, this variation of angular po- 
sitions brings them face to face, increasing their Coulomb 
repulsion in the axial direction. The comparison with 
simulations run in the pseudo-potential (equation [2|) 
confirms that this behaviour is induced by the rf motion. 
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FIG. 5. (color online) Distorsion of the double-ring structure 
viewed in the (z, j/)-plane (left picture), and from a point of 
view rotated of — 7r/4 around the i/-axis. The structure is 
stationary over the rf-period time scale. This picture is taken 
when the distance between the two centers of mass is about 
2 fim. 



The evolution of the distortion during the increase of 
LOz explains the increase of standard deviation cr^ roughly 
situated between 0.75 MHz and 0.93 MHz (Figure [J). 
When the transition frequency is reached, the cm. of 
the two structures match roughly, but the standard de- 
viation takes its maximal value, showing that the struc- 
ture is not yet a planar ring. For higher values of w^, 
the standard deviation shows a steep downward slope, 
until a frequency value slightly lower than 1 MHz. The 
structure can then be considered as a planar single-ring 
structure. 

The distortion of the structure described above and 
the coupling between the tangential and axial compo- 
nent of motion it induces, explain the initial discrepancy 
between the radial and tangential temperatures when the 
structure is a double-ring (see figure [3]). The axial mo- 
tion being not subjected to micro-motion, it can reach 
lower temperature than the motions in the radial plane. 
The coupling between axial and tangential components 
of the motion allows to dissipate a part of the tangential 
thermal kinetic energy through the axial cooling. As the 
frequency lUz is increased, the temperature increases 
progressively and meets the temperature associated to 
the radial directions, T^. At the same time, the axial 
temperature Tz shows a steeper drop between 1 MHz 
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and 1.25 MHz, reaching the Doppler hmit temperature 
at the frequency where the two other temperatures join. 
This behaviour suggests that the axial component of the 
motion of the ions is then decoupled from the two trans- 
verse directions, like expected for a full two dimensional 
system [lol |. 

The comparison between Figures S] and [3] shows that 
for this small 100-ion system, informations on the dimen- 
sionality of the system given by the variation of the ions' 
positions, cr^, and the temperatures of the system do not 
coincide. The decoupling observed on the temperature 
arises when the parameter already shows a planar 
structure. As a consequence, each of these parameters 
may be a relevant criteria to determine the structural 
phase transition. 

IV. ABOVE THE ADIABATIC LIMIT 

The single ring structure presented in the preced- 
ing section seems very promising for applications such 
high precision spectroscopy, quantum information or for 
metrological applications [9|, [l^l . For all of these, it is 
mandatory to reduce the amplitude of motion of the ions 
to below the chosen transition wavelength (the Lamb- 
Dicke regime [l^) to cancel the first order Doppler effect. 
For trapped ions, Doppler laser-cooling is sufficient if the 
trapping frequency along the laser propagation direction 
is large enough. In practice, for calcium ions cooled to 
the Doppler limit, a motional frequency of 1 MHz allows 
the ions to reach the Lamb-Dicke regime for an optical 
transition. One of the main drawbacks of the ring struc- 
ture compared to its linear counterpart, the ion chain, 
found in quadrupole traps is the non null rf-driven motion 
(micro-motion) at their equilibrium positions. Figure [3] 
shows that in the single ring configuration, the motion 
along the axis is decoupled from the motion in the ra- 
dial plane where micro-motion, and thus rf-heating, take 
place. The Lamb-Dicke regime can thus be reached for 
laser propagating along the trap axis. 

To quantify the impact of the micro-motion on the 
ion dynamics and laser cooling, the Mathieu parameter 
qx^ relevant for this issue in quadrupole traps, can be 
extrapolated to multipole trap and is then related to the 
more general concept of adiabaticity of the trajectories 
[l[ . The adiabaticity parameter iqad depends on the radial 
position of the ions inside the trap and at the potential 
minimum i?min, it takes the simple expression [6|: 

Vad (i?min) = 2 Vfc - 1 W./f! . (8) 

It has been demonstrated, numerically and experimen- 
tally in a 22-pole 0, [H [lol, that for -qad (r) < 0.36 the 
trajectories are stable but few studies concern this adi- 
abaticity parameter in the context of cold samples. In 
this section, we focus on the role of the micro-motion 
on the self-organization of the ions and relate it to the 
adiabaticity parameter. Figure [H] shows the single ring 
structure formed by 100 ions located at the same radial 



potential minimum and trapped by an axial potential 
characterized by w^/27r equal to 1 MHz (rjad = 0.173), 
2 MHz {j]ad = 0.346) and 3 MHz {j]ad = 0.520). For 
increasing a;^, we observe a bunching of the ions at the 
inter-electrode angular positions. 
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FIG. 6. View in the radial plane of structures obtained for an 
equilibrium radius i?min = 120 and, from left to right, ujz = 
1 MHz, 2 MHz and 3 MHz. From top to bottom, the pictures 
correspond to three phase of the rf: cos (fit) = —1; 0; 1. 

This bunching is not observed in the pseudo-potential 
approach of the trapping electric field and is attributed to 
the micro-motion of the ions. Some insight in the role of 
micro-motion can be gained by considering a cold sample 
sitting close to i?min, where the radial trapping potential 
can be approximated by a local harmonic potential. As- 
suming a null potential for Raim it can be written like 

f/ioc (i?min + Sr) = ^m {Srf , (9) 

with the characteristic frequency Wioc = \/k — 2 uj^ @- 
The ion motion is assumed to be due only to the rf-driven 
micro-motion whose amplitude can be written like: 

£,,{Rmin,'P,t) = (10) 

(-Rmin) II (cos (4(p) f - sin (Aip) (p) cos {Qt) , 

where the amplitude of the micro-motion He*^ (i?min) || is 
easily deduced from equation [T] Only the radial contri- 
bution of this motion engenders a significative potential 
energy variation during the rf-period. Its time average 
depends on the angular position (p like 

t^(i?min,<p) UjI^^ lie;, (i?Mn) |P cos^ (4^) (11) 

and could explain the bunching observed on figure [51 To 
test this assumption, we have created a configuration 



where only two ions bunch in one of the local angular 
minimum potential localized around tpo = {'^^ + l)7r/8. 
To that purpose, we have simulated a set of 9 ions in 
the rf-potential of the same octupole trap, with ujz/2n 
initially set at 0.5 MHz, and slowly increased to 4 MHz, 
keeping i?min = 120 /im constant by increasing the rf po- 
tential. The final stable configuration is then made of 
seven ions and a unique pair of ions, localized at each 
intcr-clcctrodc angular positions. 

Assuming a local description ~muj^g^S/2 for the force 
responsible for the bunching, the measure of the rela- 
tive distance S of the ions forming the pair allows one to 
deduce the characteristic frequency uj^es from the equi- 
librium condition with the Coulomb repulsion: 
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(12) 
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One can also evaluate the force by deriving the time- 
averaged micro-motion induced potential Uioc {Rmin,^), 
for an angle (p ^ ipo + Sip : 



Rm 



Sip (p (13) 



Writing like —muj'^Rai^Sip, one can connect the char- 
acteristic frequency to the local radial potential. The 
results of numerical simulations show that the predicted 
bunching frequency uj^ is larger than the measured fre- 
quency oJmes by a factor ^/2, all along the simulation in- 
creasing ujz at constant i?min- The disagreement by a 
factor y/2 is also observed when the same simulations are 
run with 13 ions in a dodecapole trap. According to the 
results obtained from the simulations, the local bunching 
caracteristic frequency is given by: 



{Rmin) II 
Rjnin 



(14) 



Eq. ll4l shows how the bunching frequency depends on the 
characteristic trapping parameter and ring configuration. 
Indeed, using the dependence of the micro-motion with 
the position in the trap, one can show that uJip scales like 
iol/V,. Hence, the bunching force scales like R^i^uj^ / fl'^ . 
The dependence of this bunching force with uJz was il- 
lustrated on figure [6] for fixed R^in- Its dependence with 
-Rmin for fixed uJz can be seen on figure [7] and confirms 
an increase of the bunching with Rmin- These different 
studies show that the so-called adiabaticity parameter 
r]ad is not the relevant criteria to estimate the role of 
micro-motion on the ion organization. Indeed, the four 
structures plotted on figure [7] share the same rjad- This 
behaviour shows that in the general description of the ion 
dynamics in a multipole trap, micro-motion has a non- 
negligible influence not only on the thermal propoerties 
of the sample but also on its structural properties, and 
has thus to be taken into account. 



FIG. 7. View in the radial plane of a 100-ion structure, for 
ujz/2n — 2 MHz, and R^in ~ 120 fim, 140 ^m, 160 /im and 
180 ^m, varied by modifying the rf-potential amplitude. The 
picture corresponds to a null rf-phase. 



V. CONCLUSION 

In this paper, we have numerically studied the double- 
ring to single-ring structural transition of a Doppler- 
cooled 100-ion system in the rf potential of an octupole 
trap. The simulations showed that, in the double-ring 
configuration, the temperatures associated to each com- 
ponent of the motion expressed in cylindrical coordinates 
are different. In our initial configuration, the tangen- 
tial temperature is lower than the radial one by almost 
one order of magnitude. This is attributed to a coupling 
of the tangential motion to the axial one, which is ef- 
ficiently Doppler-cooled due to the absence of rf-driven 
motion along the axial direction. When the double-ring 
structure is compressed, the cm. of each ring merge 
at the predicted transition frequency but the structure 
shows a distorsion which is maximal at this frequency. 
This behaviour is attributed to the micro-motion of the 
ions, inducing an axial Coulomb repulsion which can be 
considered as a pseudo-potential between the ions of the 
different rings. A planar single-ring structure is thus ob- 
tained for a higher value of the predicted transition fre- 
quency. When UJz is further increased, the axial tempera- 
ture reaches the Doppler limit value showing a decoupling 
of the axial and transverse components of the motion. 

For a cold ion system, we have also shown that it is 
possible to increase the adiabaticity parameter beyond 
the empirical limit, keeping stable ion confinement. At 
a fixed ring radius, increasing the adiabaticity parame- 
ter results in a bunching of the structure, for which an 
expression of the local bunching frequency is proposed. 
We have also shown that changing the ring radius at 
constant value of the adiabaticity parameter leads to the 
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emergence of the bunching phenomenon. 
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